Dynamical exponents of an even-parity-conserving contact process with diffusion 
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We provide finite-size scaling estimates for tfie dynamical critical exponent of tlie even parity- 
conserving universality class of critical behavior through exact numerical diagonalizations of the 
time evolution operator of an even-parity-conserving contact process. Our data seem to indicate 
that upon the introduction of a small diffusion rate in the process its critical behavior crosses 
over to that of the directed percolation universality class. A brief discussion of the many-sector 
decomposability of parity-conserving contact processes is presented. 
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I. INTRODUCTION 

One of the most important theoretical programs in 
modern statistical mechanics is that of the understand- 
ing and classification of the critical behavior of nonequi- 
librium interacting particle systems. A first major step 
towards this program was given by the so-called directed 
percolation (DP) conjecture ^^, which in its original 
formulation stated that all dynamically driven continu- 
ous phase transitions about a single absorbing state in 
single-component systems with a scalar order parame- 
ter in the absence of internal symmetries are in the DP 
universality class of critical behavior [pi . In this form the 
conjecture has been confirmed in a host of model systems, 
among others the basic contact process M, Schlogl's 
models for autocatalytic chemical reactions P;pl, and a 
phenomenological Euclidean field theory of high energy 
hadronic collision processes (Reggeon field theory) pM . 
Further investigation revealed that the DP universal- 
ity class accommodates moregeneral models, with more 
than a single component ||8|-[l0|| as well as with multiple, 
in some cases infinitely many absorbing states (O-nsl. 
Actually, even some non-equilibrium models without ab- 
sorbing states at all were found to share some of the DP 
exponents p6|"pJ[. 

It has been found, however, that not all phase transi- 
tions involving an absorbing state fall into the DP univer- 
sality class. Examples of models exhibiting non-DP crit- 
ical behavior range from probabilistic cellular automata 
pO| and two-temperature kinetic Ising models [|T|, to 
interacting monomer-dimer 123] and monomer- monomer 
models P3| , branching and annihilating random walks 
with even number of offspring |2^-p7[|, and a class of 
parity-conserving contact processes that is the subject 
of this paper 1 28 1 . The common feature of all these mod- 
els is that the number of interacting particles, whether 
they are actually particles or are domain walls, is locally 
conserved modulo 2, i.e., their local dynamical rules con- 
serve parity. It thus appeared at first that local conser- 



vation laws were affecting the critical behavior of non- 
equilibrium systems, as one may have guessed from his 
knowledge of equilibrium critical phenomena. The new 
universality class that emerged became known as the par- 
ity conserving (PC) universality class. Controversy arose 
when some parity-conserving models were shown to be- 
long to the DP universality class. In [2^], it was noticed 
that when one adds a parity-conserving external field to 
the dynamics of a certain monomer-dimer model with two 
equivalent absorbing states for which the number of do- 
main walls is conserved modulo 2, the universality class 
of the model crosses from the PC class to the DP class. 
This kind of behavior was then subsequently observed in 
other models [^,pl|, and nowadays it is believed that 
it is not the symmetry of parity conservation that deter- 
mines the critical behavior in these systems, but that it is 
the presence of two equivalent, Z2-symmetric absorbing 
states that matters. This has led some people to refer to 
the new universality class as the directed Ising universal- 
ity class, in allusion to the fact that the Ising model has 
two Z2-symmetric equivalent ground states. Whether the 
PC class is as robust as the DP class remains largely an 
open question, and there is lot of room left to research; 
for a review, see I32U33]. 

In this paper we provide finite-size scaling estimates 
for the dynamical critical exponent z = h'\\/i'± of the 
PC universality class of critical behavior through exact 
numerical diagonalizations of the master operator of a 
certain parity-conserving contact process in one dimen- 
sion [^. We show that upon the introduction of dif- 
fusion in the model its critical behavior crosses over to 
another universality class that seem to be characterized 
by DP exponents, although our data are not very good 
in this case. We also present a brief discussion of the 
many-sector decomposability of the state space of a class 
of parity-conserving contact processes, and how diffusion 
breaks this structure, making the models with diffusion 
strong candidates to show new types of dynamical critical 
behavior. 



II. PARITY-CONSERVING CONTACT 
PROCESSES WITH DIFFUSION 

Contact processes may be viewed as models for the 
spread of an epidemic among individuals living in a d- 
dimensional lattice. We consider a parity-conserving gen- 
eralization of the basic contact process first introduced 
in Ig^. In this class of processes, called CP(m), an ar- 
ray of m adjacent healthy individuals surrounded by 
k infected individuals X becomes infected at rate fcA, 
while an array of m adjacent infected individuals becomes 
healthy at unit rate. Pictorially, in one dimension we 

have for CP(1) the elementary processes X00 -^ XX^, 

mX A 9XX, X%X ^ XXX, and X i 0, correspond- 
ing to the usual basic contact process, while for CP(2) 

XXX^, 

1 n.A 



we have the elementary processes X000 

000X A %XXX, xmx ^ XXXX, and XX 

and analogously for m > 2. We clearly see that CP(to) 

processes conserve the number of particles modulo m. 

For finite systems with an absorbing state, the steady 
state coincides with it. For the CP(1) process, it is simply 
given by the completely empty lattice = (0,0, .. . ,0). 
CP(2) has, besides 0, the two other absorbing states 
(0,X,0,X,...,0,X) and (X,0,X, 0, . . . ,X,0). Under 
periodic boundary conditions these two configurations 
are the same; let us call it 0. Whether belongs to 
the space of allowed configurations depends on the par- 
ity of the lattice size L and on the parity sector TV mod 2, 
where N is the number of particles of the initial config- 
uration. If L is odd, then the only absorbing state for 
CP(2) is 0, since it is impossible to fill the lattice with 
a repetition of 0X's only. For a lattice of even size L, 
the state has L/2 particles, and this number has to 
be compatible with the parity of N mod 2: if L = 4/, 
with / e N, and N mod 2 = 0, then is an absorbing 
state of the process; if otherwise L ^ 41 + 2, then is an 
absorbing state only if TV mod 2 = 1. In this work, for 
definiteness we carried out our calculations on lattices of 
odd size L and on the N mod 2 = parity sector, such 
that the only steady state of the process is 0. Our meth- 
ods do of course equally apply to the other sectors of the 
dynamics as well. 

The existence of more than one absorbing state for 
CP(2) already signals a general feature of CP(m) pro- 
cesses, that the number of absorbing states of these mod- 
els grows with m. Actually, for m > 3 the number 
of absorbing states grows exponentially with the sys- 
tem size L, a property called many-sector decompos- 
ability first verified in a class of adsorption-desorption 
processes of m-mers on the lattice p3,p^. The number 
Im{L) of absorbing, fully jammed states in these sys- 
tems grows asymptotically as Im[L) ^ 2(j)^ , where (p is 
the largest real root of (p™ = 20™^^ — 1; for m — 3, 
4> = \{1 -f \/5) ~ 1.618, the golden mean, while for 
TO = 4, ~ 1.839. Notice that 4> <2 for all to < oo. This 
many-sector structure of the phase space can, however, 



be broken by adding diffusion X% ^ %X at a finite rate 
/i to the processes, since it allows the otherwise jammed 
states to evolve, eventually leading to a state with an 
array of to, adjacent healthy or infected sites that can 
then react according to the CP(to.) rules. Diffusion thus 
reduces the number of absorbing states of CP(m) from 
Im{L) ^ 20^ to Im{L) = m, corresponding to the num- 
ber of TO-parity equivalent absorbing states. While it is 
well known that diffusion, as long as the diffusion con- 
stant remains finite, is an irrelevant perturbation for the 
basic contact process CP(1) [p6i, for systems with more 
than one absorbing state diffusion may become a highly 
relevant perturbation, changing the critical behavior of 
the process, as recentl y v erified in the pair contact pro- 
cess with diffusion [BTlBSl. 



III. FINITE-SIZE SCALING 



As is well known [|39|-|41|, we may write the master 
equation for reaction-diffusion processes on the lattice as 
a Schrodinger-like equation in Euclidean time. 



dt 



\Pit)) = ~H\Pit)), 



(1) 



with \P{t)) the generating vector of the probabilities 
P{n,t) = {n\P(t)) of observing the particular configu- 
ration n = (ni,n2, . . . jUl) G {0, X}^ at instant t, and 
with H the infinitesimal generator of the Markov semi- 



(1) p(0) 



p(i) 



group. The lowest gap E\^' —E\^' = E\^' va. the spectrum 
of H may be used to perform a finite-size scaling analy- 
sis in the same way as one does in equilibrium problems 
[g|,|2[. Around the critical point A > A*, the correlation 
lengths of the infinite system behave like 



C|i cx^lcx(A-A*)-'^« cx(A-A*)- 



(2) 



where ^|| and ^^ are the correlation lengths respectively 
in the time and space directions, i/|| and v±^ are the cor- 
responding critical exponents, and z = vn /v±^ is the dy- 
namical critical exponent. For finite systems of size L, 
we expect that 



e|^i=L-^-$(|A-Al|LV-^.-), 



(3) 



where zl and v±.l are the finite versions of z and v±, 
and <I>(m) is a scaling function with $(u 3> 1) '^ u"" . 
On general grounds one expects linii^ooPi, ^l, ^j_,l = 
p* , z, u^. From Eqs. (|2[) and (^) we obtain 

HkA^DIkMK)] _ ln[C||.L"(AI)/C||,L(AI)] 



\n{L/U) 



\n{L"/L) 



Zl, 



(4) 



which through the comparison of three different system 
sizes L' < L < L" furnishes simultaneously A^ and zl- 
Of course, ^||^ and the gap Ef^ of H are related by 
e,7i=Re{i?«}. 



IV. DYNAMICAL CRITICAL EXPONENTS 

We calculated the gaps of H through the power 
method, which requires only matrix-by-vector multipli- 
cations that can be carried out efficiently, and does not 
require a diagonalization in the usual, 'QR' sense, a step 
that may lessen the quality of the data. The particu- 
lar implementation of the power method we use takes 
full advantage of the presence of absorbing states in the 
process, and is also suitable for the investigation of time 
dependent properties of finite-state Markov chains |43| . 

Our results for A* and z for the CP(2) both without 
diffusion and with symmetric diffusion X% ^ %X at rate 
/I = 0.05 are summarized in Tabic |. The choice fi = 0.05 
is arbitrary, except that it represents only a small per- 
turbation to the main coupling A near the critical points, 
but not so small as to render too much unbalanced, 'stiff' 
matrices that may become difficult to diagonalize. We 
were able to diagonalize systems with up to L = 25 sites 
in a reasonable amount of time and memory space, com- 
parable to recent density-matrix renormalization group 
studies of similar systems |37| . The extrapolated values 
for the /i = case were obtained through a Bulirsch- 
Stoer extrapolation [Q, while those for the fJ. ^ case 
were obtained through a least-squares fit to the curve 
Vl — ao + aiL~^ + a2L~^ , and can be seen in Fig. |l|. The 
uncertainties associated with our extrapolated numbers 
are mainly due to finite-size effects and corrections to 
scaling, as well as to the extrapolation procedures itself. 
We obtained in the diffusionless case A* — 0.88 ± 0.01 
{lubst = 0.924) and z = 1.75 ±0.01 (cjbst = 1.802), and 
for the diffusive case A* = 0.25 ± 0.05 and z = 1.3 ± 0.1. 
(A linear ffi y^ = ao + aiL~^ to the last four points of the 
data with diffusion furnishes A* = 0.38 ±0.04 with a cor- 
relation coefficient 7 = 0.989, and z — 1.36 ± 0.06 with 
a correlation coefficient 7 = 0.989.) We thus see that 
the pure CP(2) without diffusion has a critical behav- 
ior governed by the dynamical critical exponent of the 
parity-conserving universality class, for which the best 
known value to date is given by zpc = 1-750 ± 0.005) 
p5[ . The value of the critical point agrees well with the 
value A* = 0.8935 ±0.0004 found in ^. Upon the intro- 
duction of diffusion, however, the number of steady states 
of CP(2) reduces from two to one, and according to our 
data its critical behavior crosses over to that of another 
universality class. In a first moment the critical exponent 
seem to be converging to z ~ 2. Data for larger lattice 
sizes, however, point toward a lower value of z, definitely 
different from ^; = 2, at least as far as our finite-size data 
go. We also observed a non-monotonic behavior of the 
data, that unfortunately seems to be more common in 
this type of calculation than it would be desirable. This 
non-monotonic behavior is probably related with the ex- 
istence of a whole critical (in the present case, diffusion- 
like) phase to the right of the critical point, and has 
already been observed in other studies of systems with 
extended critical phases |lq;p7| . Our extrapolation gives 



TABLE L Finite-size data and extrapolated values for the 
critical point and the dynamical critical exponent z of CP(2) 
with diffusion X0 ^ 0X at rate fi. The numbers between 
parentheses represent the estimated errors in the last digit of 
the data. 



System sizes 
L',L,L" 


AI 


= 

Zl 


AI 


0.05 

ZL 


7,9,11 


0.584 505 


1.916 330 


0.705 564 


1.975 081 


9,11,13 


0.607 946 


1.898 353 


0.759 294 


2.013 974 


11,13,15 
13,15,17 


0.625 558 
0.641431 


1.884 819 
1.872 620 


0.791 092 
0.798 318 


2.046 900 
2.056 315 


15,17,19 


0.656 023 


1.861476 


0.783 969 


2.034 662 


17,19,21 
19,21,23 


0.669 350 
0.681 446 


1.851398 
1.842 353 


0.753 899 
0.715 444 


1.985 477 
1.921 154 


21,23,25 


0.692 390 


1.834 263 


0.675 905 


1.856 800 


Extrapolated 


0.88(1) 


1.75(1) 


0.25(5) 


1.3(1) 



z = 1.3 ± 0.1, closer to the exponent of the DP univer- 
sality class, namely z^p ~ 1.5807 |[43| , [l5[| , than to other 
known values. A DP critical behavior is what one could 
have expected on the basis of the DP conjecture, since the 
process with diffusion has a single absorbing state. Intu- 
ition, however, has proved not to be reliable in guessing 
the critical behavior of nonequilibrium systems with lo- 
cal conservation laws or many absorbing states, as one 
infers from the recent history in the field, and the pres- 
ence of an additional symmetry might well have driven 
the critical exponents of CP(2) with diffusion to those of 
the PC universality class. Diffusion is then seen to be a 
highly relevant perturbation for systems with more than 
one absorbing state. 



V. SUMMARY AND CONCLUSIONS 

In summary, we have conducted a finite-size scaling 
analysis of the dynamical critical exponent of an even- 
parity-conserving contact process through exact numer- 
ical diagonalizations of its time evolution operator. We 
showed that the critical behavior of the 'pure' model 
with two absorbing configurations belongs to the parity- 
conserving universality class. In the presence of symmet- 
ric diffusion, however, the number of absorbing config- 
urations in the model reduces from two to one, and its 
critical behavior crosses over to that of another universal- 
ity class that seem to be that of the directed percolation 
process. This behavior is in accordance with what one 
expects on the basis of the directed percolation conjec- 
ture, at least for small diffusion rates. 

It would be interesting to perform time-dependent 
Monte Carlo simulations of the CP(2) model both with 
and without diffusion in order to determine its critical 
exponents more precisely, as well as to investigate the 
critical behavior of other members of the CP(r7i) class of 
processes. In particular, it would be very interesting to 
investigate the CP(3), since it was argued on the basis of 
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FIG. 1. Finite-size data and least-squares fits for tlie crit- 
ical point A* and tlie dynamical critical exponent z of tlie 
CP(2) model with diffusion. 



Monte Carlo simulations and Pade approximants analy- 
ses that this process does not suffer an ordinary second 
order phase transition [g^. Moreover, the CP(3) pro- 
cess has the many-sector decomposability property, and 
it would be interesting to see how diffusion changes the 
scenario for the phase transitions, if any, in this model. 

It seems that the question as to the roles of symmetries 
and many-sector decomposability on the critical behavior 
of nonequilibrium interacting particle systems is far from 
being answered, and more numerical and analytical work 
has to be done before a consistent scenario emerges. 
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